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O ' Abstract 
(N . 

The STIT tessellation process was introduced and examined by Mecke, 
^ I Nagel and Weifi; many of its main characteristics are contained in [5]. 

In [3], Mecke introduced another process in discrete time. With a geometric 
distribution whose parameter depends on the time, he reaches a continuous- 
time model. In his Conjecture 3, he assumed this continuous-time model to 
' be equivalent to STIT. 

Ph . In the present paper, that conjecture is proven. An interesting relation arises 

[ to a continuous-time version of the equally-likely model classified by Cowan in 

' p|. This will also clarify how Mecke's model works as a process in continuous 

time. 

MSG (2000): 60D05 
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Mecke, Nagel and WeiB developed a model for a tessellation process in a window W 
^ '. in M.'^,d > 2, which has the characteristic property of being stable unter iteration. 

In [5], many of the characteristics of the so called STIT process were examined and 
shown. 

In Mecke introduced a process in discrete time in which at every discrete time 
■ step a decision is taken whether to change the state of the tessellation process or 

^ . not. This decision depends on the number of quasi-cells (i.e. cells of the tessellation 

and other "virtual" units that have no interior and cannot be hit by a line but can 
nonetheless be chosen for "division") in a given tessellation and on whether a line 
which is distributed independently of the chosen quasi-cell hits the cell or not. The 
tessellation process changes its state only if a ("real") cell is chosen for division and 
the line hits that cell and thus splits it in two cells. Otherwise, the line is said to be 
dismissed or rejected. 

Mecke then introduced a random variable which is geometrically distributed with a 
parameter that depends on a time t G [0, oo). He then related this discrete random 
variable with his discrete-time process. This yields a continuous-time model for 
which it is not immediately clear how it works as a process. In his Conjecture 3 in 
[1], he considers that his continuous-time model should be equivalent to the STIT 
model in continuous time. It is the main goal of the present paper to prove this 
conjecture. 
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To this end, in section [21 the STIT tessellation in continuous time is introduced 
in a way that is convenient for the further considerations. Some definitions are 
given which are used throughout the present paper and a key property of the STIT 
tessellation, the distribution of the n-th jump time given a certain cell configuration 
sequence, is examined. 

In section [3l the Mecke process in discrete time is introduced formally and certain 
properties of this process, e.g. the discrete waiting time for a state change and the 
distribution of the i-th jump time of the process, are examined and proven given 
the existence of a certain cell configuration sequence. 

In section m Mecke's way of turning his discrete-time process into a continuous-time 
model is examined closely for a given cell configuration sequence. The cornerstone 
of this section is the proof of the identity of distributions (in STIT and Mecke's 
continuous-time model) of the number of jumps until a given time t G [0, oo) under 
a given cell configuration sequence. 

While in section Othe identity of cell configuration distributions and thus the identity 
of the tessellation distributions for a certain time t G [0,cxd) of STIT and Mecke's 
continuous-time model will be shown it turns out that there is no obvious way of 
understanding Mecke's continuous-time model as a process. 

This problem is addressed in section El Here, one of the models Cowan introduced 
in [2] is examined. It turns out that the number of jumps in Cowan's equally-likely 
model has the same distribution Mecke employs for the number of decisions in his 
transition from discrete to continuous time. It turns out that one can plug the 
process of the jumps in Cowan's model into Mecke's discrete-time process to get 
Mecke's continuous-time process which by then has been shown to be equivalent to 
the STIT process in continuous time. 

2 The STIT tessellation 

Let in this paper the window W G M."^ always be a convex and compact polygon in 
the plane with a non-empty interior. 

2.1 Intuitive introduction 

The intuitive idea of the STIT construction introduced in [3 Section 2.1] allows the 
following notation of a STIT model. Throughout this paper, the planar case will be 
examined. 

Let ["Hji^] be the measurable space of all lines in where the a-algebra is induced 
by the Borel cr- algebra on a parameter space of "H. For a set A C we define 

[A] = {ge'H:gnAj^iD}. 

Let A be a non-zero locally finite and translation-invariant measure on ["Hji^] which 
is not concentrated on one direction. For A, < A([W^]) < oo is true. 
It should be noted that we will be working with 'lifetimes' of cells that are condi- 
tioned on the existence of exactly these cells. As the probability for the existence 
of a certain cell composition is zero, one has to retreat to regular versions of condi- 
tioned probabilities which are rigorously defined in e.g. The intuitive handling 
we make use of here remains correct. 

The initial window shall have a lifetime ~ £^(A([VF])), i.e. its lifetime is ex- 
ponentially distributed with a parameter A([VF]). At the end of its lifetime, the 
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window is divided into two cells by the segment within W ol a line 71 whose distri- 
bution is A{[W])-^A{- n [W]). 

Let us call the resulting cells Ci and C2. Let us further have exponentially dis- 
tributed independent random variables Ti,T2, ... with Tj ~ £^(1), j = 1,2, .... Then 

the cells Ci and C2 shall have lifetimes = j^j^^ti ~ £(A([Ci])) and = 
^(^([^2]))- At the end of each lifetime, the segment within Cj of a line 
j£, distributed A([Cj])~^A(- fl [Cj]), divides the cell whose lifetime is over, resulting 
in another two cells, and so on. 

So at every point in time, the cell Cj from the extant cells Ci, has a lifetime 

■^'^j ^ A{[c-]) '^^ ^ ^(^([Q])) with some independent of all other Tj. 
One of the immediate consequences of this is that for the waiting time for the 
state of the whole tessellation in W to change (i.e. the time span until one of its 
current cells is split) 

T,^£(A([Ci]) + ... + A([Cfc])) (1) 

holds. 

We denote as Y^{t, W) the state of the tessellation in the window W at the time t. 
If, at the time t, there are cells Ci, Cm in the tessellation, we define the random 
closed set 

m 

Y^{t,W) = \JdCj\dW 

where dCj denotes the boundary of a set Cj. 

Then, the STIT process is given by W):t>0). 

2.2 Preliminary considerations 

It the following, many probabilities will be conditioned on the existence of a certain 
cell configuration. 

Definition 1 Let Cn he the set of cell configurations from the time t — ( with W 
being the only cell) until before the i-th jump time in a tessellation process: 

Ce = {{M^}, {Ci,i, (^1,2}, {C2,i, C2,2, (^2,3}, {Ce-1,1, Ce-i^e}}, 

while the order of the cells Ci k after the i-th jump time is irrelevant. It holds 
Ci = {{Co,i}} = {{W}}. 

A sequence {Ce : i & N) is called compatible if (for all £ E N) there are indices 
ki, ke^2 (ind k[, k'^_2 such that Ce-2,kj = Q-i,*;^. for j = 1, — 2, and for the 
remaining indices ke-i,k[_^ and k'^ the equation 

holds. 

Let further be for k — 1,2, ... 

It will turn out to be easier to work with each C(, being a set (instead of an £-tuple) ; 
this does not constitute a problem for ordering because, for every k — 1,2, 
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there is only one set in Ce, that has exactly k components. 

Obviously, 1 < Ck < and Ck < k almost surely for A; > 1. In the following, all 
sequences (C^ : G N) are considered compatible. 

Let us now look at a tessellation process in the window W and a certain cell config- 
uration sequence (C^ : £ G N) on which we condition. At the time t = 0, the window 
is empty. It takes a time span (waiting time) Ti (depending on W) until a first 
segment is thrown into the window. After another waiting time T2 (depending on 
the cells after the first jump) which is independent of all other waiting times (under 
the given condition), another segment is thrown into the window, and so on. The 
£-th jump time is defined by = Yl^j=i where the waiting times Tj are all 
independent of each other under the given condition. If at a time t there are £ cells 
Ce-i^i, Ce-i^i then we again denote 



Y{t, W) = \J dCe-i,, \ dW 
j=i 

and {Y{t, W) : t > 0) as the process. 

If we have a cell configuration sequence (Q : £ G N) but examine a property that 
depends only on the first n cell configurations (such as the n-th jump time t„), it is 
sufficient to consider the cell configuration C„. 

Under the assumption of the compatibility of the (Q : £ G N), the following lemma 
holds. 

Lemma 1 For the n-th jump time t^ of the STIT tessellation process in W under 
the condition of the cell configuration Cn and for n = 1, 2, ... the equation 

n ^ 

p(tf < t|c„) = 1 + (-1)" 5^ e-^mw -Q (2) 



k=\ iG{l,...,n}\{fc} 



holds. 



Proof 

The lemma will be proven by induction. 

The base case for n = 1 is correct because 

P(tf < t\{{W}}) = P(tf <t) = l- e-^([^l)* = 1 + (_i)ie-^([^])^i* ■ 1 
and £1 = 1. 

Equation (|2]) is equivalent to the probabihty density ft„ (conditioned on Cn) being 
written as 

ftSx) = f P(t„ < x\Cn) = (-1)"+^ 5^ A([W^])Ae-^([^l)^^^ n T^.- 

^ k=l ie{l,...,n}\{A:} ^ ' 

Let us use the abbreviation 



n^"'= n c._c 

ie{l,...,n}\W 



Then 

n 

fM = (-i)"+^$^A([t^])Anr)e-^(['^i)^'=^ 



k=l 
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is equivalent to ([2]). It will be shown that stepping from tn to tn+i via the waiting 
time Tn+i and the convolution formula yields the same result as equation ([2]) applied 
for tn+i directly. 

According to the convolution formula (which is applicable because of the conditional 
independence of all waiting times of each other), for the density ft„+^{x) 

/t„+i (x) = / ftn («)/r„+i (x - u)du 
Jo 

is correct. Because of Tn+i ~ £{A{[W])Cn+i) = SiJ2]j=i ^ii^j])) per equation 
([I]), for the density /t„+i 

/r„+,(x -u) = A([iy])£„+ie-^([^])^"+^(^-") 

holds. So it is sufficient to show that 

= ^^+1 A([W^])£fcni"+')e~^([^])^^- 

(3) 

holds. Let us begin with a consideration arising from interpolation theory: 

Let us have the function f{x) = 1 for all a: G M. Then, this function / = 1 can be 

interpolated as the linear combination p of Lagrange polynomials (see also [Ij) 

n n 

pw=e/(-') n n 

k=\ ie{l,..,n}\{fc} " * fc=l iG{l,..,n}\{fc} " * 

with being arbitrary interpolation points. This function / = 1 is a polyno- 

mial itself, having degree which is less than the number of interpolation points ^— 1, 
thus the linear combination of Lagrange polynomials p is f itself: f{x) = p{x) = 1 
for all X G M. Let us now denote Xk = Ck, with the Ck as in the lemma. Then we 
get the equation 

k=l ie{l,..,n}\{k} ^ * 

or in the special case x = Cn+i, 



k=l le{l,..,ri}\{fe} ^'^ 

From this, we get the following sequence of equivalent equations. 
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En TT -Cn+i-£i 
k=l LLie{l,..,n}\{k} Ck-Ci 

Z^fc=l Cn+i^Ck llie{l,..,n}\{fc} Ck~C, 

Z^k=l Ck-Cn+1 llie{l,..,n}\{fc} Ck-C, 

„-A([w])Cn+ix V" r, n^"+^) 



n" i 

_n" i 

iii=l £„+i-£i 
/- TT(n+l) 

-'-n+l-L-'-ri+l e 
_^n+l ^ T-r(n+l) -A([H^])£;,x 

2^k=l ^k^^k ^ 
_ spn+1 ^ TT(n+l) A([H/])£fcX 



^ l^k=l ^ki-i-k l^k=l '-k^^k ^ 

g-A([w-])£„+ix Y2-1 Ckllf^^\l - e'^([^l)('C"+i-'Cfc)^) 

iV^(/l/JJ^.„+ie Z^A:=l'^fc^^fc A([iy]){£fc-£„+i) 

Thus, equation ([3]) is true. □ 

3 The Mecke process in discrete time 
3.1 Introduction 

In [3], Mecke develops a new process in discrete time. 

There are lines 7j , J = 1, 2, that are i.i.d. according to the law A([M^])^^A(-n[W^]). 
Further let us use, independently of 7j, independent aj, j = 1,2, ... where aj is uni- 
formly distributed on the set {1, 

If a line 7 does not contain the origin o then 7+ shall be the open halfplane 
bounded by 7 which contains the origin. Correspondingly, 7" is the open half- 
plane bounded by 7 which does not contain the origin. As the distribution of 7 is 
assumed translation- invariant, we can neglect the possibility of 7 going through the 
origin as the probability of this is zero. 

Let be Co,! = W, Ci^i = Wn-ft and Ci,2 = Wn 7f . Forn = 2, 3, ... we define 

Cn-i,j if j e {l,...,n},j an 

Cn,j = < Cn-l,an H 7n if 3 = "n 

. Cn-i,a„ n 7+ if j = n + 1 

These entities Cnj are called quasi-cells. Some of these quasi-cells are empty. 
Those quasi-cells that are not empty will be called cells. 

From this, we can deduce a random process: After each decision time n, n = 
1,2,..., we consider the tessellation Tn consisting of the quasi-cells Cn,ij ■■■^Cn^n+i- 
This decision time is called the n-th decision time accordingly. If, at that decision 
time, the number of cells (i.e. non-empty quasi-cells) actually changes, that decision 
time is called a jump time. Obviously, the k-th jump time is that decision time at 
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which the number of cells reaches k+1. Let us denote the random closed set of the 
closure of the union of cell boundaries that are not part of the window's boundary 
at a step n for the tessellation Tn as 



n+l 

Yf(n,W)^[jdCr.,\dW. 

Then {Yj^{n, W) : n e N) is called the Mecke process in discrete time. 

3.2 Discrete waiting time for a state to change 

In the following, we will denote cells as Cj, i.e. without the tilde which is used when 
examining quasi-cells. 

Lemma 2 After n — 1 time steps, let k cells Ci, C^. be extant in the tessellation in 
the window W (n = 1,2, k = 1,2, ...,n). The probability that at the n-th decision 
time a fixed cell Cj e {Ci, Ck} is split is 

fiAt the n-th decision time the cell C,- is splitlCi, Ck) = ttttttIt"- 

A([M/J) n 

The probability that at the n-th decision time any of the cells C\, ...,Ck is split can 
be calculated as 



' A([Q]) 



¥{At the n-th decision time any of the cells Ci, Ck is split\Ci, Ck) — — \ {mn\ 

j^i ^\[^\) 

Proof 

The probability that a line 7„ hits the cell Cj is x^j^jy- Additionally, independent 
of the line, one of the n quasi-cells will be selected for division. Only if the line 7„ 
hits the selected cell, there is a change in the state of the tessellation process. So, 

P(At the n-th decision time a cell is split |Ci, Ck) 
— X]^=iF(At the n-th decision time the cell Cj is split |Ci, Ck) 

_ A([C,]) 1 

Z^j=l A{[W]) n 

_ is^k A([Cj]) 
n Z^j=l A{[W]) ■ 

□ 

If, after the (n — l)-th decision time, we have k cells, i.e. non-empty quasi-cells, in 
the tessellation Tn-i, this state will be denoted Y(^n-i),k- Let us denote -'^y(„_i) ^ as 
the discrete waiting time for this state to change or the discrete jump time 
span between the (A; — l)-th and the A;-th jump time, i.e. the time span until 
a line segment actually falls into a cell. Let us further denote the k cells Ci, Ck- 

Lemma 3 For i — 1,2, n — 2,3, ... and 2 < k < n the equation 

(n-1)! T{n + e-Ck-l) 



{n + i-l)\ r(n-A) 
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holds. For n = 1, 

nxvo,, = 1) - 1. 

The measure defined in this way is, for all k and n with 1 < k < n, a probability 
measure onN\ {0}. 

Proof 

Atn— 1, the window W is hit by the hne almost surely; therefore its hfetime (which 
here equals the waiting time) is almost surely exactly one time step. 
Let n be greater than 1. Then, the waiting time for a state K„_i,fc to change is i 
if and only if the tessellation is not changed in the n-th, the (n + l)-th until the 
{n + i — 2)-th time step but eventually jumps into another state in the {n + i— l)-th 
time step. The probability for the tessallation to be unchanged in the {n + j)-th 
time step is 1 — ^^Ck, the probability for it to jump into another state in the 
{n + £ — l)-th time step is n+e-i - follows 

P(Xy(„_,,,=£|Ci,...,C,) 
^ n+i-l^k 11^=0 (-'- ~ ^^k) 
n+i-l'-'k i-i-3=0 n+j 

_ 1 p r{n+l-l-Ck) 
~ n+i-l'-'k [n+i-2)\ T(n-Ck) 

_ p (n-1)! r(n+£-£fc-l) 
- >~'k{nU-\)\ r(n-£fc) • 
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Afterwards, it follows with Fubini's theorem for all n G N \ {0, 1} and for all k < n 
P(l<^y(„_„,, <oo|Ci,...,C,) 

= E,~iP(^y(„-,),.=^|Ci,...,C,) 

Eoo ^ (n-1)! r{n+e-Ck-l) 
£=1 '^fc(„+^_l)! r(n-£fc) 

Eoo r{n+l-Ck-l) 
e=i {n+e-iy. 



(c) 



(n- 


-IV 




-jCk) 


u ^ 


-1)! 




-jCk) 


(n- 


-1)! 






in- 


-1)! 




-jCk) 




-1)! 


•'=r(n 






1 


'kr{n 






1 


■'=r(n 






1 



(4) 



(71—1)1 /-uw ^-tj.—L.i.— \ _t / 1 1 in,!) \ 



In (a), r(n) = J^^-^e-^dt. In (b), 

x"-^e-^da; = (n - l)!e-* ^ - 

k=0 

for n e N is the incomplete gamma function. In (c), partial integration is used and 
the fact ^ 

j- (r(n) - r(n, t))^^ f x^'-^e-'^dx = T-^e"*. 
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In (d), 

limt^ot-''^ (r(n)-r(n,t)) 



limt^^o 



= 0, 

due to I'Hopital's rule and n > almost surely because of n > 1 here. □ 
Special cases include 

P(^y(„_i),fe = 1|Ci,...,Ca:) = -Ck 

and 



V{Xy^„_,^^^ = 2\C^, Ck) = (1 - Uk):^Ck. 



3.3 The discrete jump times Xg 

Let us denote with Xi, i = 1,2,..., the discrete jump times, i.e. the time span 
from the beginning of the observation in the window W at the time until the £-th 
change of the state of the tessellation process. 

For the discrete jump times Xi and X2 the following distributions hold: 

P(Xi = 1) = 1 
and for n = 2, 3, ... it follows from Lemma |3] 

P(X2 = n|C2) = P(Xy,,, = n- l|Ci,i, = ^2^ r(2-£2) ' 

Due to the Markov property of the Mecke process of tessellations in W and the fact 
that the shape and size of the cells 6*2,1, 6*2,2 and 6*2,3 does not depend on the jump 
time X2, it follows for n = 3,4, ... 

n-l 

P(X3 = nICa) = J2 P(^2 = fc|C2)P(Xy, 3 =n- k\C2,i, 6*2,2, 6*2,3) 

k=2 

with 

63 = 62 U {{62,1, 62,2, 6*2,3}} = {{W}, {61,1, 6*1,2}, {^2,1, 62,2, 62,3}}. 

In general, for i = 1, 2, ... and n = i, i + 1, ... we get 

n-l 

P(X, = n\Ce) = J2 ^(^^-1 = ^l^i-iM^YKi = ^ - k\Ce-i,i, 6,_i,,). (6) 
k=e-i 
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Because of the compatibility of the sequence (C^ : £ G N), 



Note that while neither shape nor size of the cells depend on the jump time, this is 
not true vice versa because the jump times do depend on the shape and size of the 
cells. 

We can derive 

Lemma 4 For the Mecke process in discrete time and for i = 2, 3, n = i, £+1, 



P(X, = n\Ct) = 
holds. This is a probability measure, i.e. 

oo 



I T{n-Ci) -j-r 1 



, T(2-Ci) J- J- A 

t=2 \ ^ je{2,...,e}\{i} 1 



(7) 



for all i = 2,3, ... and all cell configuration sequences {Ci : i E N) . 
Proof 

The proof is by induction. 

For the base case {£ = 2), we have 



n\ ^r(2-£2) n\ ^1(2 -A) 

which is equivalent to equation (It follows from convention that HieO^* ~ 
In the induction step it is sufficient to show that under the assumption of correctness 
of the lemma for £ and inserting it into equation ([6]) on one hand and inserting £ + 1 
into equation ([7]) on the other hand, both equations yield the same result. Let us 
use the abbreviation 

n zhr- 

ie{2,.../}\{i} ' ^ 

and begin with our consideration arising from interpolation theory again: 
Let us have a function / : M — )■ M with £ — 1 values f{x2), f{xi) at interpolation 
points X2, ■.■,xi. Then, this function / can be interpolated by the linear combination 
p of Lagrange polynomials (see also [I]) 



P{.x) = ^f{xi) JJ 



i=2 je{2,..,n\{i} ^ 

If the function / is a polynomial itself, having a degree which is less than the number 
of interpolation points i — 1, the linear combination of Lagrange polynomials p is f 
itself: f{x) = p{x) for all x G M. Let us consider the function 

^ ^ ^ k=2 k=2 
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This function is a polynomial of degree i — 2. Thus, we can say 



f{x)=v{x)=Y,j{x,) n 



i=2 jG{2,../}\{i} ^ 

Let us now denote = Ck, with the as in the lemma. Then we get the equation 

^ ^ i=2 j€{2,..,i}\{i} ^ 



or in the special case x = C 



^ «=2 je{2,..,e}\{i} J 

Multiplying equation ^ by ^Yit-xT+i) ^^"^ having /(x) = p|^~^j as above, we get 

Er(2 — r(£ — Ci) T-r — _ 

Separating the fraction in the product and multiplying by y^^tc^i g^t 

^ ^ r(2 — Cf^j^i) r(£ — c ) I T— I- 1 I 1 ( T~r 



Dividing both sides by (^11^=2 ~ ^dji using the abbreviation ^ and multi- 
plying by r(2-£^+i j ' equation becomes 

r(?2-£<;+i) r(£-£j) 2,^+1 ^ r(n - £^+1) 2/+! 
r(2-£,) r(£-£,+i) ^ r(2-£,+i) ' 

Subtracting a term on both sides yields 

sr^l V{n-Ci+x) T-r2,^+l Tji-Ci) _ s^l T(n-Ci) T-r2,£+l 
l^i=2 V{2-Ci) llj r(€-£i+i) l^i=2 V{2-Ci) Hi 



r(n-£<,+i) T-r2,£+l _ V^^ / r(n-£0 T-r2/+l 
■r(2-£f+i) Vu+i l^i=2 yr{2-Ci) Hi 



v^£+l / r(n-A) T-r2/+l 
Z^i=2 I r(2-£i) Hi 
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Thus, we get 



_ sr^e+l ( r(n-£,) T-r2/+l\ 

2^i=2 I r(2-A) lli ) 



S-^l r(n-£f+i) T-r2/ T(i-Ci) sr^i 1 T-r2,£ T(n-Ci) 

2^i=2 r(2-£,) Hi (A-£f+i)r(^-£,+i) Z^i=2 r(2-£,) Hi 



^ r(n-£^+i) Yr2/ T(l~C^) v-^ V{n-Ci+r) T-r2,£ T{n-Ci) 

2^1=2 r(2-£,) Hi (£,-£,+i)r(^-£,+i) Z^i=2 r(2-A) Hi (£,-£<+i)r(n-/:i+i) 



l^i=2 

(a) 



E 



i=2 



r(2-£,) 



(£,-£f+i)r(£-£,+i) (A-£^+i)r(n-£^+i) y_ 

r{2-Ci) Hi l2^fc=£ r(fc+i-£f+i) 



En— 1 
fc=£ 2^ 



i=2 \r{2-Ci) ili 



rik+i-Ci+i) ■ 



Multiplying both ends of this chain of equations by f Y[i=2 '^0 ' ^^t 



v^^ ( r(fc-£,) T-r2/ 
^i=2 \r{2-Ci) Hi 



€+1 



r(n-£f+i) 

r(fc+i-£f+i) 



/"r-ri'+l ^ [v^+l f r(n-£,) T-r2,^+l 
\^lli=2'^*y [Z^i=2 \^r(2-£i) 11 



Finally, by multiplying by (-l)'^^g, we get 



k\ r(n-£f+i) 
n! r(fc+l-£f+i) 



= (-ir'ir(nS£.)[ES(?g5gnr' 

This, however, is equivalent to equation ([7]) 



1 / r(n-£0 T-r2,^+l 



n- 



i=2 \^r(2-£i) Ili 

and thus the induction step is finished. 

The equation (a) above is shown itself by induction; it will be shown that 



n-l 

E 

k=t 



r(n- A) 



(10) 
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From this, (a) follows immediately. 

The base case n — 1 = £ is straightforward: 



1 r(^~£,) i-c, 

Vlyt-Ci) 1 



r(£+i-i;,+i) • 

Let equation f llUp be correct for n — 1; it follows for n 

En r{k—Ci) 
k=i r{k+l-Ce+i) 

^ fs^n~l r(k^C^) \ I r(n-£,) 

y2^k=e rik+i-Ci+i) J r(n+i-£,+i) 



C,-Ce+i \r{e-Ce+i) r{n-Ce+i) J ' r{n+l-Ci+i) 

rji-CQ r(«-A) . {Ci-i:e+i)r{n-i:i) 



1 / r{e-C^) _ r{n-Ci)in-Ce+i-C^+Ce+i)) 



1 / r{e-d) _ r{n-Ci){n-Ci)) 



rje-CQ rjn+i-Ci) 



This is the result equation fITU]) yields directly for n. Thus the equation is true, also 
for every i due to possible index shift. 

To show that F{X£ = -ICe) is a probability measure on {£,£+ 1,...}, let us first 
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consider that similar to equation 



Eoo 1 r(n~£,;) _ 1 Y^oo _!_ r°° j.n-Ci-l„-tJj. 

n=l n\ r(2-£i) ~ V{2-Ci) l^n=l n\ JO ^ ^ 



r(2-A) 



1 r°° ^-t+~Ci-i„t 



TOO _f, 

r(2~£~) Jo ^ V"^ 



e* 1 



r(^,t) 



r(2-£,)rw Jo 



1 

r(2-£0rw 



r(2-£i)rw ^ io ^' *'^'^) 
r(2-£i)rw tJ^^^ ~ 

T{t-x) 



holds. Then, with f{x) = p|^_^j as above 



2 \ r{2-c,) iiie{2,...,£}\{j} C,-Cj 
oo 1 r(n— £i) 



^1)^ (ni=2 Z]j=2 (nje{2,.../}\{j} £~~C7 ) ^n=e ^ r(2-£i) 



-ir(nl2A) EI2 (n 



iG{2,...,n\{i} ; £ir(£)r(2-£i) 



r(^-£.) 



-'-J (i-iy. 2^i=2 yilje{2,...,e}\{i} c,-c, J r(2-£,) 



■^Yji^i-^y-'EL (n 



j€{2,...,e}\{i} c, -Cj J ^{2-c^) 



1 r(<;-o) 



-1)! r(2-o) 



(£-1)! 
= 1. 

Thus, the lemma is proven. 



□ 



4 The Mecke model in continuous time 
4.1 The distribution of ^{t) 

In [U Section 4], Mecke introduces a mixed line tessellation model such that the 
tessellation T* at the continuous time t E [0, 00) corresponds to the tessellation Tu{t) 
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at the discrete random time z/(t) where for the distribution of z/(t) 

F{u{t) = A;) = e-* (l - e"*)*' ,k^O,l, ... 
holds. For general A([H^]), the distribution is 

F{u{t) = fc) = e-^([^l)* (1 - e-^([^l)*)' ,k^O,l, ... 

This is the geometric distribution with parameter e"^^'^'^*. A possible interpretation 
here is that the decision times are no longer at equidistant discrete times n — 1, 2, ... 
Instead, the law describes how many decisions take place until the time t. The 
i>{t) are assumed independent of all other random variables that are used in the 
construction of the Mecke process. 

Corollary 1 For the probability that until the time t there have taken place at least 
n decisions, the equation 

niy{t) > n) = (1 - e-^([^l)*)" 

holds for n — 1,2, ... 
Proof 

oo oo 

niy{t) > n) ^J2^{iy{t) = i) = ^e-^([^])* (l - e-^^I^^)*)' = (l - e-^^I^D*)" . 

i=n i=n 

□ 

Let us now denote the number of jumps until the time t in the Mecke model as 
ri^{t), i.e. the number of jumps until the discrete time ^{t) in the discrete-time 
Mecke process. 

Corollary 2 For the conditional probability that, under the condition of a cell con- 
figuration Ct, until the time t, there have taken place at least i jumps, for i — 1, 2, ... 

oo 

P(r/^(t) > e\Ce) = 5] (1 - e-^([^l)*)" P(X, = n\Ce) (11) 
n=e 

holds. 

Note that due to F{r]^\t) > i\Ce) = 1 - ¥{r]^{t) < i\Ce) one does not need any 
knowledge of the jumps after the i-th. 



4.2 Comparison with the properties of STIT 

We will now compare the properties of the continuous-time Mecke model with those 
of the STIT process. Generally, properties relating to the STIT process arc denoted 
by the upper index ^ , properties relating to the Mecke model are denoted by the 
upper index ^ . 

Obviously, for STIT the equation 

holds with tf the £-th jump time and ri^{t) the number of jumps until the time t. 
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Lemma 5 Under the condition of the cell configuration Ce, for £ = 1,2, ... 

holds for all t G [0, oo). 
Proof 

For e = l,Ci = {{W}} holds. With respect to STIT, 

P(tf < = 1 - e"^([^])* 
according to Lemma [H For Mecke's model examined here, 

P(r7^(t) > l|Ci) = (1 - e-^([^])*)¥(X, = l|Ci) = (1 - e-'^^t^l)*) -1 = 1- e"^"^!)* 

holds because of P(Xi = 1) = 1. 
For ^ = 2, 3, ... it is to show that 



E~ , (1 - e-^'W'T (-1)'^ (nl.A) [El. (S^S n 



1 



2 \^r(2-£0 iiie{2,...,^}\{i} 



(12) 

In this equation, the left-hand side is won by inserting the equation from Lemma H] 
into equation (ITTi) ; the right-hand side follows from Lemma [H 
The left-hand side of equation (fT2l) can be re-written as 

- I ^) l^i=2'~'^\iij<^{2,...I}\{i} C^-Cj ) 2^n=t\^ ^ ) n\T{2-C,) 



which remains to be shown. With A = 1 — e ^(['^D* one gets 

1 r{n-Ci) 
n\ r(2-£i) 

^-A(fVKl)t^" 

n! r(2-£,) 



^ _ (>-M[w])ty j_r(n-£i) 

Y^oo An 1 r(n-£i) _ ^^£-1 _ A([VK])A" 1 r(n-£0 
Z^n=0^ n! r(2-£i) Z^n=0 l"^ ^ - 



r(2-£i) Jo Z^n=0 n\ ""^ Z^n=0 ^"^ J n! r(2-£i) 

- 1 r°° p-a;^-A-lpAx J -AfWlHN" 1 r(n-£,) 

r(2-£i) Jo ^ e ux Z^„=o ^ n! r(2-£i) 

r(2-£i) Jo ^ "-^ Z^n=0 y n! r{2-Ci) 

substitution 1 r°o „-« <^ u 1 ^ _ _ A([VK])A" 1 rjn-Ci) 

r{2-Ci) Jo ^ \\-AJ I-A"" l^n=0\^ ^ ] n\T(2-Hi) 

~ r(2--£i)^ ^ Z^n=0 l.-^ n! r(2-£i) 

~ r(2-£i)^ ^ ^ ^n=0V ^ J n! r(2-£i) 

_ e-A([w])£,t _ , _ ^~A([T^])A'^ 1 r(n-£,) 

(£i-l)£i Z^n=0 J n!r(2-£i)- 
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Re- writing this as 



in the right-hand side of the original equation f[T^ one gets the new equation 



c-n^+iv^ r /^n A y^^-i fi ^-Affi^pA'^ i r(«-£,) 

i + (_i)V^(m)*ni2T^ 



(13) 



which is equivalent to the original equation ( IT2|) . We will now show that his equation 
is true. 

In (IT^ . because of 



fc=0 ^ ' 



.;^^fcg-A([v^^])fc^ 



there appear the summands e ^(['^l)^*^/^ = 0,...,£ — 1. With respect to the sum- 
mands, we will compare the coefficients. These can be written as 



^ ,fn\ r(n-£.) _ fe/A {i-h)Y{l-L,) 
> \k)v{2-C,)n\ ^ > VV^!(^-A)r(2-A) 



(14) 



for A; = 0, 1, 1 and z = 2, 3, 

First, the correctness of equation (|T4|) is shown by induction, keeping k and z fixed. 
It is quite obvious that it is equivalent to show (IT^ or the simplified equation 

^ V{n - £,) _ 1 r(f - A) 

^ (n-fc)! (£-A;-l)! ■ ^ ^ 

n=k 

The base case, for £ = 2, is correct for the two possible values of A;, namely A; = 
and k = 1: For /c = 0, we get 

r(-A) + r(i-A) = - ^^^~^^^^ +r(i-A) = -lr(i-£,)(-i + A) = ^^^Zcr 

For A; = 1, we get 

r(i - A) ^ 1 r(2 - A) 

1! (2-1-1)! 1 -A • 
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Now, the induction step follows. Let us assume that equation (1T5|) is correct for (. 
and any k < ^ — 1. Then, for £ + 1 we get 



r(n-£,) 
--k {n-ky. 



Z^n=k {n~ky. (f-fc)! 



1 r{e-Ci) , rje-c) 
{i-k-i)\ k-d {e-ky. 



r(^-A) / e-k , 1 

(^-fc)! I A:-£i 



r(^-£.) 



- A; + A; - A 



(£-fc)!(A:-£,) 
(^-fc)!(fc-£i)'^^ 

^ r(£+i-A) 

(^+l-fe-l)!(fc-£i) 

which is what equation (IT^ says for £ + 1. For k = i, we have 
r(£-£,) 1 r(£+l-£. 



0! (^ + 1-^-1)! 



r(£-£.). 



Thus, equations (fTSjl and consequently (|T4|) are true. 

Let us now make use of this result and examine the coefficients. 

The coefficient for A; = on the left-hand side of equation (ITS]) is 



1)^^^ J2i=2 (nje{2,...,^}\{i} ) ( (o) 

J-J Z^i=2 ^^llje{2,...,£}\{i} (^-l)!r{2-A) 



e\ i£-o)T(e-Ci) 
e\{o-c,)r{2-Ci) 



i&{2,...,e}\{i} C,-Cj ) (^-l)!r{2-£,) 

r(^-o) 



(^-l)!r(2-0) 



_ 

- (^-1)! 

= 1. 

This corresponds to the coefficient on the right-hand side of equation ( IT^ . Here 
was, once again, exploited the fact that the interpolation polynomial of a polynomial 
of a certain degree (here i — 2) that is less than the number of known data points 
(here i — turns out the be that original polynomial. The interpolation polynomial 
was evaluated at the point 0. 

In similar fashion, it is shown that the coefficients k = 2, 3, — 1 are because 
(k-c )r{2-c ) again such a polynomial of a degree {i — 2) less than the number of 
known data points {£ — 1). Here, the interpolation polynomial is evaluated at the 



19 



point as well. 



1^1=2 yi-i-jG{2,...,i}\{i} Ci~Cj J \k) mk-Ci)V{2~C^) 
l^i=2 \ iij&{2,...,i}\{i} Ci-Cj ) \k) 



kfe\ {i-k)r{e-Ci) 
e\{k-c,)r{2-Ci) 



^2e-i+ksr^ (rt o--Cj \ fe\ c^{e^k)r{e-Ci) 

e\{k-c^)r{2-Ci) 



_iy£-i+k i^e^ o(i-k)r{e-o) 



e\{k-o)r{2-o) 



0. 



This again corresponds to the coefficients of the right-hand side of equation (fT3|l . 
Finally, the coefficient for k = 1 is evaluated. If the equation 



■J 



holds, then the coefficients are equal. Beginning with the obvious 1 = 1, we get 



{e-iy. 

(e-i)T{e~i) 
(^-i)!r(2-i) 

EL (n, 

y ^) 2^i=2\nje{2,...,e}\{i} c-Cj J '^aii-Ci)r{2-c^) - 11 



e{2,.../}\{i} c,-cj (^-i)!r{2-A) 



j=2 

and thus the desired equation. Here, the interpolation polynomial was evaluated at 
the point 1. 

Thus, all coefficients of the e~^([^l)''* are equal for k = 0,1, — 1. So, the equa- 
tions are equivalent. □ 
While it follows from this result that the numbers of jumps have identical distribu- 
tions for any given t G [0, oo) and a sequence of cell configurations (C^ : £ G N) in 
Mecke's continuous-time model and STIT (in continuous time) it is not yet clear 
how the Mecke model describes a random process in continuous time in which, at 
certain points of time, the state of the tessellation process is changed. Up to now, 
for every point in time t G [0, oo), the number z/(t) of decision steps was evaluated 
separately according to 

P(zy(t) = k) = e^\l - e-'f 
or, for A([M^]) ^ 1 according to 

P(z.(t) = k) = e-^"^])* (1 - e-^([^])*)' (16) 

respectively. Then the discrete process was observed until this number v{t) was 
reached. It is not yet clear, however, how the process comes from a time ti > to 
another time t2 > ti because the random numbers ui = z/(ti) and U2 = ^^(^2) are not 
independent in such a scenario. For instance, 1^2 > J^i must hold necessarily. This 
problem will be addressed a bit later. 
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5 Equivalence of STIT and Mecke's continuous- 
time model 

Let us first generalize the result of Lemma [5] for a fixed sequence of cell configura- 
tions. 

Lemma 6 The cell configurations Cf^ of Mecke's models (in discrete or continuous 
time) and Cf of the STIT model have, for £ = 1,2, an identical distribution. 

Proof 

The proof follows from induction with respect to the algorithm according to which 
the cell configuration develops. 
Obviously, = Cf = {{W}}. 

Let now Cf^ = Cf . So, after the {(, — l)-th jump time there are exactly £ cells, equal 
in both configurations. The selection probabalities of a cell to be split in the ^-th 
step are equal; the probability for the A;-th cell to be selected is 

P(The cell C^-i,fc is selected for division |C£_i^i, ...,Ci^i/^ — ^ ^''^^ 



E}=iA([^^.-i.]) 

for Mecke's as well as for the STIT model. For the STIT model this was shown and 
used rather frequently, for Mecke's model this is true because of 



n the £-th division step the cell C^-i,fc is split |C^_i,i, Ci-i^i) 
En=i nCi-i,k n 7n ^ 0|C,_i,i, C,_i,,, U ... U n 7„ ^ 0)P(X, = n 

n=e p{(c,_i,iU...uc,_i,^)n7n7^0|c,_i,i,...,Q_i,,)^'^"^'^ - "'\^<^) 

s-^oD A([C^-i,fc]) nA{[W]) W>(X — r,\r ^ 

2^n=e nA{[W]) A([C,_i.i])+...+A{[C,_i,f])^'^"^^ " '^\^<^) 

Eoo A([C^^i,fc]) Jtsf Y — r^ln \ 
n=£ A([C,_i,i])+...+A([Q_i,,])^l"^^ - ^l'-^) 



A([C,,i.fc]) 



j:7=enXi = n\Ce) 



A([C,_i,i])+...+A([Q_i,«]) 

^ A([g^-i,fc]) 

A([Cf-i,i])+-+A([Q_i,,])- 

The final equation is true because the measure P(X^ = -jC^) is a probability measure 
as proven in Lemma HI 

The division of a cell follows the same division rule, namely according to a A law, 
so that the distribution of both configurations is the same. □ 
From this it follows immediately 

Theorem 1 Mecke 's continuous time model and STIT have identical distributions 
for a fixed t G [0, oo) . 

Proof 

For every fixed cell configuration Ce the identity of the conditional probability for 
Mecke and STIT follows from Lemma [51 The distribution of the cell configurations is 
identical for both models according to the above Lemma El Thus, for all i = 1,2, ... 

{V^{t) >i} = {v^{t) > £} 
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for all t G [0, oo). Because the cell configurations have an identical distribution for 
every i = 1,2, ... and are independent of the jump time, the distributions of the two 
models are identical for every t G [0, oo). □ 
Thus, Conjecture 3 in is proven. 

6 The Mecke process in continuous time 

At the end of section H] it was mentioned that it is not yet clear how to understand 
Mecke's continuous-time model as a process. This section is to show a solution for 
the problem; we will find a connection to one of Cowan's models. 

6.1 Cowan's equally-likely model 

Cowan summarized eight different models for cell division (in discrete time), coining 
the terms of selection rule and division rule, in [2]. One of the four selection rules 
he introduced is the equally-likely rule. In this, in a given tessellation with n extant 
cells, the probability of a certain cell to be selected for division is ^, without regard 
of its perimeter or area or anything else. Cowan mentions a rather straightforward 
way of extending this model towards continuous time. 

A way to do this is to give each cell a lifetime which is exponentially distributed 
with a fixed parameter, say 1, and independent of all other lifetimes. At the end of 
each cell's lifetime, this cell is divided according to a division rule (here the division 
rule according to the law A([C])~^A(- fl [C]) for a cell C is suitable), resulting in two 
cells each with a lifetime exponentially distributed with the given parameter 1 and 
independent of all other lifetimes. 

Thus, when there are k cells in a given tessellation (i.e. after k — 1 divisions), 
the lifetime of each cell is £^(l)-distributed, and the whole state of the tessellation 
process has a waiting time for change that is £^(A;)-distributed, as mentioned above. 
Let be the n-th jump time and Nf the number of jumps up until time t. Let 
Tf ~ £{k) for A; = 1, 2, .... Then 

k=l 

and 

n 

= max{n : E rf < t}. 

k=l 

{N^ : t > 0) is denoted as the process of the number of jumps in Cowan's equally- 
likely model. 

Let us first take a look at the two following general lemmas. 

Lemma 7 Let n G N \ {0} be fixed. Let further Sn be the sum of independent 
exponentially distributed random variables Ti, T„ with Tj ~ S{j) for j = 1,2, ...,n. 
Then ^ 

P(5'„ <t)= [ ne-"^(e^ - 1)"-Mx = e-"*(e* - 1)" 
Jo 

holds. 
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Proof 

Again, the proof is by induction. For n = 1, obviously 

F{Si <t)= P(Ti <t)= [ 1 ■ e"^-^ ■ (e^ - ifdx = [ e-^'dx = 1-6"* 

Jo Jo 

holds which is true according to the condition Ti ~ £^(1). 

Let the lemma be true for n. Then, because of Sn+i = Sn+Tn+i with T„+i ~ £{n+l) 
and the independence of 5*^ and Tn+i, for the density of Sn+i 

fs„+^ (x) = .fs^+T^+^ {x) = fs„ {u)fT„+, (x - u)du 

= (n + l)e-("+^)^ ne"(e" - 

= (n + l)e-("+^)^[(e" - l)"]:i=S 

= (n + l)e"("+i)^(e^ - l)*" 
holds. Integration delivers the second equation in the lemma. □ 



Lemma 8 Let Nt = max{T2, : Yl^=i — ^} denote the number of jumps after 
waiting times Tj, j = 1, ...,n, until the time t. Then for k = 0,1,2, ... 

F{Nt = k) = e-' {I - e-'f 

holds. 

Proof From the distribution of the Sk, k = 1,2, one gets 

P(iVt = k) 

= ¥{Sk<t<Sk+i) 

= nSk<t)-F{Sk+i<t) 

= e-^* (e* - 1)^ - e-*^*e-* (e* - 1)*= (e* - 1) 

= e-'=*(e*-l)''(l-e-*(e*-l)) 

= e-^^ (e* - 1)'' (1 - 1 + e"*) 

= e-('=+i)* (e* - if 

= e-*(l-e-*)\ 

For A'^t = 0, the result follows from Lemma [7] immediately. □ 
It is straightforward to relate these two lemmas with the corresponding properties 
of Cowan's equally-likely process: The n-th jump time t^ in Cowan's process corre- 
sponds to Sn from Lemma [71 the number of jumps Nf until the time t in Cowan's 
process corresponds to Nf from Lemma [HI 
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6.2 Relation between the models 

If, in the process of the Cowan equally-hkely jump times, the parameter is not 1 but 
A ([14^]) one gets 

P(Xi = A:) = e-^([^])*(l-e-^(['^l)*)' 

for the equation from Lemma [HI This, however, is exactly equation (|T6l) which is 
the generahzation of the equation Mecke employs for the number of decision steps 
until the time t where Xt in Mecke's notation is called u and in this paper usually 

Thus, Mecke's model can be understood in such a way that at each time t the waiting 
time for the state to change until a new quasi-cell, not necessarily a new cell arises 
is a random variable with exponential distribution whose parameter is the number 
of quasi-cells at that time, multiplied by the factor A([VF]). 

Put differently: If at the time t there exist n quasi-cells in a Mecke model (the 
number of real cells is irrelevant) then this state (of n quasi-cells) has a pseudo- 
waiting time f;^ ^ £{nA{[W])). 

Generally, the number of quasi-cells corresponds to the number of decisions plus one 
(in Mecke's process); the number of cells corresponds to the number of jumps plus 
one (in the process of Cowan equally-likely jump times). 

Lemma 9 The number of decisions in the Mecke model with time t G [0,oo) has 
the same distribution as the number of jumps in Cowan's equally-likely model if the 
parameter A{[W]) is equal: 

X? ^ u{t). 

The distributions of the pseudo -waiting time T*^ of the state between the {n — l)-th 
and the n-th decision time in Mecke 's process and of the waiting time T^ of the state 
between the {n — l)-th and the n-th cell division in Cowan's equally-likely model are 
identical. It holds 

for each n = 1, 2, ... 

Definition 2 Let us have a window W CR"^. Let {Yj^^{n, W) : n eN) be the Mecke 
process in discrete time as described in section\^ Let {N^ : t > 0) be the process of 
the number of jumps in Cowan's equally-likely model. Then for every t G [0, oo) we 
define 

y,^(t,M^) = r/(A^f,t^) 

and the Mecke process in continuous time as (y/^(t, W) : t > 0). 

Finally, due to the Markov property, the identity of the one-dimensional distributions 
(Theorem [H Lemma [9]) and the identity of the stochastic kernels as per Lemma El 
we have the following 

Theorem 2 Let us have a window W C Let {Y^{t,W) : t > 0) be the STIT 
process in continuous time as described in section IE o-nd {Y^'^ {t,W) : t > 0) the 
Mecke process in continuous time as described in definition Then 

(y^(t, W):t>{)) = {Y^\t, W):t>{)). 
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